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A new method is proposed for a treatment of ideal quantum gases in the microcanonical ensemble 
near the thermodynamic limit. The method allows rigorous asymptotic calculations of the average 
number of particles and particle number fluctuations in the microcanonical ensemble. It gives 
also the finite-volume corrections due to exact energy conservation for the total average number 
of particles and for higher moments of the particle number distribution in a system approaching 
V") the thermodynamic limit. A present consideration confirms our previous findings that the scaled 

variance for particle number fluctuations in the microcanonical ensemble is different from that in 
' the grand canonical ensemble even in the thermodynamic limit. 
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C3 ' I. INTRODUCTION 

The statistical hadron gas approach to nucleus- nucleus (A+A) collisions (see e.g. Ref. [l( and recent review 0) is 
■ rather successful in describing the data on particle multiplicities in a wide range of the collision energies. Usually one 
\ considers a thermal system created in A+A collision in the grand canonical ensemble (GCE) . The canonical ensemble 
t-H ■ (CE) 0, 0, IS El 0, El or even the microcanonical ensemble (MCE) has been used in order to describe the pp, 
pp and e + e~ collisions when a small number of secondary particles are produced. In all these cases, the statistical 
ON . systems are far away from the thermodynamic limit, so the statistical ensembles are not equivalent and exact charge 
conservation or both energy and charge conservation laws have to be taken into account. The CE is relevant also for 
systems with a large number of produced particles, e.g., a large number of resultant pions or large nucleon number 
in p+A collisions, but a small (of the order of 1 or smaller) number of carriers of conserved charges, such as strange 
hadrons [jj, antibaryons Q or charmed hadrons Q. This may happen not only in elementary but also in p+A or 
f^i . even A+A collisions. 

""p^ ' The analysis of the fluctuations is a useful tool to study the properties of the system created during high energy 
+3 , particle and nuclear collisions (see, e.g., Refs. 0, 0> EH El EH EH El E3 an d the review papers fig). In A+A 
1 1 ' collisions one prefers to use the GCE because it is the most convenient one from the technical point of view and due 
O \ to the fact that both the CE and MCE are equivalent to the GCE in the thermodynamic limit when the size of the 
2 ' system tends to infinity. However, the thermodynamic equivalence of ensembles means only that the average values 
of physical quantities calculated in different ensembles are equal to each other in the thermodynamic limit. It was 
demonstrated for the first time in Ref. |19| that multiplicity fluctuations are different in the CE and GCE even in the 
thermodynamic limit. These results have been then verified and extended in Refs. pel I2TL |22 | . The particle number 
fluctuations in the MCE have been considered in our paper [2^ and they have been shown to be different from the 
' GCE results even in the thermodynamic limit. 

In this paper we present a new method for the study of particle number distribution in the microcanonical ensemble. 
The method is based on the analysis of the asymptotic behaviour of moments of the particle number distribution. In 
contrast to the previously used microscopic correlator approach |15|. the new method is more rigorous, consistent and 
mathematically justified. It elucidates some subtleties. Along with fluctuations, it allows to calculate the finite- volume 
corrections to the thermodynamic quantities. 

The paper is organized as follows. In Section II we review the method of microscopic correlator and demonstrate 
its limitations. A detailed description of our new method is given in Section III. We summarize our consideration in 
Section IV. 
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II. MICROSCOPIC CORRELATOR 



Let consider the quantum system of non-interacting neutral particles. We review here the application of the method 
of Ref. |l5j| to the calculation of the particle number fluctuations in the systems with the exact conservation of energy 
imposed. 



2 



The GCE partition function for a single quantum state with momentum p has the form 

z p = ^cxp (-^ nj , (1) 



where T is the system temperature, e p = \/p 2 + m? and m is the particle mass. The sum in Eq. Q runs over the 
number of particles n — 0, 1 for the Fermi statistics and n = 0, 1, 2, . . . oo for the Bose statistics. Summing up the two 
terms for the Fermi statistics, or the infinite geometric series for the Bose statistics, one gets: 



(2) 



where 7 = +1 and 7 = — 1 for Bose and Fermi statistics, respectively. The GCE average values are calculated as 
(fc = l,2): 



in",), . = ^2 n k w p (n), (3) 



1p)g.c.e. 



where 



Wp (n) = Zp 1 exp( - n) (4) 



is the probability to observe n particles in the given quantum state. It is easy to see that 

M _ (~ T ) k dkz P 



in , 



For k = 1 we get the familiar Fermi (Bose) distribution 



exp 1^1 



(5) 



(0) 



and for k = 2 



2 p)g.c.e. = eXP (T) + 1 = (n p ) s . c . e .[l + (l + 7 ){n p ) 9 . c . e .] . (7) 



From Eqs. (101711 it follows: 

((An p ) 2 ) g . c . e . = (n 2 p ) g . c . e . - (n p ) 2 g c ^ = (n p ) g . c . e . [1 - j(n p ) g . c . e .} = v 2 p . (8) 

It is easy to see that 7 = in Eqs. (161811 corresponds to the Boltzmann approximation. 

Expressions © and JSJ) are microscopic in a sense that they describe the average values and fluctuation of a single 
mode with momentum p. However, the fluctuations of macroscopic quantities of the system can be determined through 
the fluctuations of these single modes. To be more precise, we will demonstrate that the fluctuations can be written 
in terms of the microscopic correlator (An p Ank) g.c.e.- This correlator can be presented as: 

(An p An k ) g . c . e . = v 2 S pk . (9) 
The variance ((AN) 2 ) = (N 2 ) - (N) 2 of the total number of particles, N = J2 P n v, equals to: 

((AjV) 2 ) 9 . c . e = ^2(n p n k ) g. c . e . - (n p )(n k ) g. c . e . = ^(An p An k ) g . c . e . = ^v 2 . (10) 

P ,k P ,k p 

We have assumed above that the quantum p-levels are non-degenerate. In fact each this level should be further 
specified by the projection of a particle spin. Thus each p-level splits into g = 2j + 1 sub-levels. It will be assumed 
that the p-summation includes all sub-levels too. This does not change the above formulation because of statistical 
independence of these quantum sub-levels. The degeneracy factor enters explicitly when one substitutes, in the 
thermodynamic limit, the summation over discrete levels by the integration: 



£...*£/*,.... (U) 
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FIG. 1: The scaled variance i l'Iti for particle number fluctuation in the grand canonical ensemble for different types of statistics. 
The lower and upper solid lines correspond to the Fermi and Bose ideal gas, respectively. The dashed line is the Boltzmann 
approximation. 



The scaled variance u in the thermodynamic limit V — > oo reads: 
((A7V 2 )) 9 . c . e . £ pfc (An p An fc ) s . c .e. 



V P 



Jo P d P V P 



Jg.c.e. 



(N) 



g.c.e. 



Jo P 2 dp (% 



The Eq. I|12l) corresponds to the particle number fluctuations in the GCE. From Eq. 1)1 2H one finds 



classical Boltzmann limit (7 = 0) 



The effects of quantum statistics lead to u> 



Jg.c.e. 



(12) 



= 1 in the 



g.c.e. 



g.c.e. < 1 for the Fermi gas (7 = —1). The strongest effect corresponds to m/T 

tt 2 /6C(3) ~ 1.368 , uj F Q Zr = T 2 /9C(3) = 



> 1 for the Bose gas (7=1) and 
0, 



0.912 



(13) 



g.c.e. I W ^V"/ ? ""g.c.e. 

and it decreases with increasing m/T (see Fig.HJ. 

The formula for the microscopic correlator is modified if we impose the exact conservation law on our equilibrated 
system. For this purpose we introduce the equilibrium probability distribution W(n p ) of the occupation numbers 
{rip}. As a first step we assume that each n p fluctuates independently according to the Gauss distribution law with 
mean square deviation v p : 



W(n p ) cx Yl exp 



(An P Y 



(14) 



To justify this assumption (see Ref. one can consider the sum of n p in small momentum volume (Ap) 3 with the 
center at p. At fixed (Ap) 3 and V — > 00 the average number of particles inside (Ap) 3 becomes large. Each particle 
configuration inside (Ap) 3 consists of (Ap) 3 • V/(2ir) 3 » 1 statistically independent terms, each with average value 
(rip) © and scaled variance v p (jHJ- From the central limit theorem it then follows that the probability distribution 
for the fluctuations inside (Ap) 3 should be Gaussian. In fact, we always convolve n p with some smooth function of p, 
so instead of writing the Gaussian distribution for the sum of n p in (Ap) 3 we can use it directly for n p . 

Now we want to impose the exact conservation laws. The conserved quantity A (the energy and/or conserved 
charge) can be written in the form A = E P a (p) n p- An exact conservation law means the restriction on the set of 
occupation numbers {n p }: only those which satisfy the condition AA = J2 P a (p)An p = can be realized. Let us 



consider exact energy conservation. Then 
because of the energy conservation as: 



A — > E (i.e. a(p) — ► e p ) and the distribution (|14|l will be modified 



W(n p ) cx Y[ ex P 



(An P 



e p An p 



dX 



IT cxp 



(An p Y 



+ iX e p An p 



(15) 
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where S (e p An p ) is the Dirac's S— function. It is convenient to generalize distribution (|15fl using further the integration 
along imaginary axis in A-space. After completing squares one gets: 



W(n p ,X) oc Yl 

exp 



(An p -\v%e p ) A 2 2 2 
2«2 + 2 Vp€p 



(16) 



and the average values (i.e. the MCE averages) are now calculated as: 

( . = ijZ^IZJk^v -W(n p ,X) 

Using Eq. (|T7|) one easily deduces 

((An, v 2 p Xe p )(An k - v 2 Xe k )) = S pk v 2 p , (A 2 ) = - «p^ . 
Therefore, one finds the MCE average for the microscopic correlator 



(17) 



((An p - ^Ae p )A) = 



(An p An k } m .c.e. = 5 pk - v 2 t p v 2 k t k (A 2 ) + (An p X) v\t k + (An k X) v 2 e p 



X 2 2 /\2\ x 2 V P € P V k €k 

pk + v p e p v k e k (A ) = 6 pk v p - — — — 

t—fp v p e p 



(18) 



By means of Eq. (|18|l one obtains: 



((AN 2 )) m . c . t 

{N) m .c.e. 



y v 2 

Z^p u p 



Ep( n p)fl-c.e. 12p( n p) g-c.e. y, p v p e p 



/ °° P 2 dp v 2 p (J °° p 2 dp v 2 e p ) 2 



Jo P 2d P 



/ °° P 2 dp (n p )g.c.e. 5™p 2 dpvle 2 p 



(19) 



The Eq. (|19|) demonstrates that the MCE fluctuations in the thermodynamic limit V — > oo can be presented in terms 
of the GCE quantities. The MCE scaled variances l|19[) for different statistics are shown as functions of m/T in Fig. [3 
The microcanonical suppression effect for the particle number fluctuations increases with the particle mass. In 



the Boltzmann approximation (7 = 0, v 2 
analytically: 



= ( n p)g.c.e. — ex p(— Ep/T 1 )) the integrals in Eq. (|T§|) can be calculated 



jf°V^ex P (-|) = Tm 2 K 2 ^) , 



2 / £p\ TTI 

p dp e p exp (^- — j = — 



p dp e p exp 
v J 



T6 



3 Wr) + Mt) - 2 ^(t) 



(20) 
(21) 
(22) 



Making use of the asymptotic behavior of the modified Hankel function K n (x) at x — > (i^o 0&) — — In a; and 
i4T„(a;) ~ 5 r(n) (§) " for n > 1) one gets in the massless limit: 



u m .c.e.(rn = 0) 



(23) 



i.e. for classical massless particles the MCE the scaled variance is quarter as large as the corresponding scaled variance 
in the GCE. For the case of Bose and Fermi statistic we obtain: 



Bose 1 n\ 



''(m = 0) = 



6£(3) 



9£(3) 



135f(3) 

2tt 4 
405£(3) 

7tt 4 



0.535 , 
0.198 . 



(24) 
(25) 
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FIG. 2: The scaled variance l|19|l for particle number fluctuation in the microcanonical ensemble for different types of statistics. 
The lower and upper solid lines correspond to the Fermi and Bose ideal gas, respectively. The dashed line is the Boltzmann 
approximation. 

The oJm.c.e. for bosons is always larger then in the Boltzmann limit and for fermions is always smaller, so we acquire 
the Bosc-cnhancement and Fermi-suppression of the fluctuations. They exist also in the GCE. However, the effects due 
to the quantum statistics in the MCE become stronger than those in the GCE (compare Eqs. (|23I25I) and Eq. (|13l) h 
As it is seen from Fig. [5] the MCE suppression effects for the fluctuations of massive particles increase with particle 
mass. One finds that with increasing m/T ratio the scaled variances for the Bose and Fermi systems approach their 

Boltzmann limit, cj m . c . e ., and from the asymptotic behavior, K n (x) ~ V^S" ex P( — ^0 + 4 " 83 T 1 ^ at x 3> 1, it follows 

that LOm.c.e. — 3/2(m/T) -2 — > at m/T — > oo. 

Comparing Eq. (|18|l and Eq. 10 one finds the changes of the microscopic correlator due to the exact energy 
conservation. First, in the MCE the fluctuations of each modep is reduced, i.e. the ((An p ) 2 ) calculated from Eq. (|18|l 
is smaller that calculated from Eq. 0. Second, the anticorrelations between different modes p ^ k absent in the 
GCE © appear in the MCE l|18l) . These changes of the microscopic correlator result in the suppression effect of the 
MCE scaled variance ui m .c.e. i|19|) in a comparison to the GCE one oj g . c . e (|12|l . 

The method considered above has several serious limitations 

• The probability distribution for An p was replaced by the Gaussian distribution 1141) . but in fact, this two 
distributions are very different. It is enough to mention that the real distribution is a discrete, and the Gaussian 
one is continuous. Although the arguments given after the equation (|14f> suggest that this replacement does not 
influence the final result in the thermodynamic limit, still, it would be nice to have a more rigorous justification 
of this fact as well as a recipe for the treatment of finite-size corrections in quantum systems. 

• The Exact energy conservation influences not only fluctuations, but also the average values (n p ) and (rip), and 
this should result in microcanonical corrections to (A) m . c . e . and (A 2 ) m . c . e . in the case of finite-volume systems. 
The approach presented in this section does not allow calculating these corrections. 

• It would be desirable to convince ourselves that the present approach gives a correct value for uj m .c.e. despite of 
the fact that it ignores the finite-volume corrections to (A) 2 ^ c e and (A 2 ) m . c . e . which are of the same order as 
fluctuations ((AA) 2 ) m . c . e . = (A 2 ) m . c . e . - (A0m.ce/ 

The rigorous method described in the next sections resolves the above issues. 
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III. THE MOMENTS OF THE PARTICLE NUMBER DISTRIBUTIONS IN THE MCE 



In this section we consider a consistent and mathematically justified method for MCE treatment of an ideal gas 
near and in the thermodynamical limit. Our aim is to find the asymptotic behavior of the thermodynamic quantities 
(the total particle number and its fluctuations) of a microcanonical thermodynamic system at large volume V in terms 
of the thermodynamic properties of a grand canonical system, which are easier to calculate. 

The methods is based on the analysis of the particle number distribution and its moments. For pedagogical purposes, 
we first apply our approach to the grand canonical system and than use the gained experience for the treatment of a 
microcanonical system near the thermodynamic limit V — > oo. 

The grand canonical partition function of the ideal quantum gas is given by a product of the grand canonical 
partition functions z p (J2J) for each quantum level p: 

Z g .o.e.(T) = Y[ Zp = IlE eX P ("^) = E Wff-ce. ({"*}) , (26) 
P P n p {n p } 

where the sum runs over all possible sets of the occupation numbers {n p }. The quantity 

W g . c . e .{{n p }) = exp (- p ^ np \ =Y[z p w p {n p ) (27) 

\ p 

is proportional to the probability to observe a given set {n p } of the occupation numbers. Here w p (n p ) is the probability 
Q to observe n p particles at the level p. 

To get the (unnormalized) probability distribution for the total particle number we multiply the above expression 

with the (5-function S (^2 p n p — N^j and sum over all n p : 



W g , c . e .(N) = E W g . c , e .({n p })6 [J2 n P ~ N ) 

{n p } \ P I 



(28) 



Now let us consider the Fourier transform (|A2(I of the above probability distribution: 

W g . c . e . {v) = exp [ivN) W g .c e . (N) . (29) 



N 

The 5-function disappears due to the summation: 



exp [ivN) 5 I rip — N ) = exp I iv n p ) , 

N \ p J \ p J 



this makes the expression (|29|) factorizable 



Wg.c.e. iy) = E exp ^E n P Wg.c.e. ({n p }) = TT Z P ^ ex P ( ivn p) w p {n p ) = \\w p (i 

{n p } \ P J p n p p 



(30) 



(31) 



Here we have introduced a Fourier transform w p {v) of the single- level probability distribution w p (n p ): 

w p (u) = exp {ivn p ) w p (n p ) . (32) 



One can rewrite the equation l|31|) as 



Wg.c.e. (f) = ex P I E l0g ^fM ) II 

V p / p 



■p ■ 



(33) 



The summation in the exponential can be replaced by integration, if the system is sufficiently large: 

gv 



Wg.c.eXv) = exp 



(2tt)3 



dTp log w p (v) 
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Let us expand the logarithm in the integral into a Taylor series: 



log w p (v) = log w p (o) + 



where, in particular, (see Appendix) 



mi(p) = (n p )g. c . e . , 

TO 2(P) = ((Hp _ ( n p)) 2 )g.c.e. = I'- 



(35) 



(36) 
(37) 



If one uses the normalized probability distribution for w p (n), then the first term in l|35|l is zero and the expression 
(|34fl takes the form 



Wg.c.e. {v) = CXp 



E 



d 3 p rrij (p) 



n 



(38) 



Now one can calculate the moments of the probability distribution W g . c . e .(N) using the values of the function 
Wg. c , e .(v) and its derivatives at v = (see Appendix). The value 



%c.e.(o)=n- 



(39) 



is used for the normalization. The first derivative is related to the average number of particles: 

1 dW g . c . e .{v) 



Similarly, 



(N 2 



iW g . c . e (0) dv 



!/=0 



1 d 2 Wg. c . e .(v) 



W g . c . e .(0) dv 2 



gV I q 

I dp(n p ) g . c . e .=N 



(2tt) 3 



drpvp 



(40) 



(41) 



The above results for the GCE are known from the textbooks. The purpose of this consideration is to demonstrate 
how our method works and to clarify our further step - a treatment of the MCE. The MCE partition function of the 
ideal quantum gas is given by 



Z m .c.e.{{n p }) = J2 S \J1 e P n P ~ E ) ■ 

{n p } \ P I 



(42) 



where the sum, similarly to Eq. I|26|) . runs over all sets of the occupation numbers {n p }. The probability to observe 
a given set {n p } is proportional to 



Wm.c.e. ({n p }) = 6 ^2 e P H P _ 



(43) 



Nothing changes if we multiply the last expression by 1 = exp(E/T) exp(—E/T): 



W m .c.e.({n p }) = exp (J^J exp (^-J^j 6 (^2 



e p n p - E . 



(44) 



The parameter T has the meaning of the temperature of a grand canonical system, which will be defined later. 
Using the properties of the (5-function, one can replace E in the second exponent by J2 P e p n p : 



W m . c . e .({n p }) = exp ( — ) exp 



■pi bp 



T 



S ^ e p n p - E 
\ p J 



(45) 
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Then the second exponential function can be rewritten using Eq. Q as 

■ y p ) =ii ex p ^n^w^ 



exp 



(46) 



where z p is the grand canonical partition function for the single p- level, w p (n p ) Q is the probability to observe n p 
particles at this level. The grand canonical system is assumed to have the same quantum levels as the microcanonical 
system under consideration. The (5-function in l)45|l can be represented as 



S I ^2 e p n p ~ E j — 2~ J ^ exp 



iX I e p n p — E 



(47) 



therefore, 



1 f E\ f + °° 

W m . c . e .({n p }) = ^ ex PyfJ J d\exp(-i\E)Y\_w p (n p )z p exp(i\n p e p ) 

/+oo ^ 
dX exp (— iXE) J^J w p (n p ) exp (iXn p e p ) , 
-oc 



where we have introduced the notation 



C = — exp 
2tt 



n 



(48) 



(49) 



The constant C is a normalization factor and is irrelevant to the physical quantities we are interested in. 

Similarly to Eq. (|28(l , we get the (unnormalized) probability distribution for the total particle number by multiplying 

the above expression with the (5-function 8 {^2 P n p ~ Nj an d summing over all possible sets {n p }: 



W m . c . e .(N) = C 



r+oo i N 

/ dXexp (—iXE) w p (n p ) exp (iXn p e p ) S I n p — N 

J °° p n p V p j 



We perform a Fourier transformation of W m . c . e . (iV): 

W m . c .e.(v) = exp (ij/iV) .Wrn.ceXiV) 



(50) 



(51) 



N 



The 5-function disappears due to (|3C)|> . and the integrand becomes factorizable: 



Wm.c.e.M = C 



C 



/+oo ^ 
(iAcxp {—iXE) JT^ w P {n p ) exp [m p (z/ + Ae p )] 
- oc _ 



P ri p 



/+oo ^ f+OO 

dAexp (— iAi?) JT^ u) p (f + Ae p ) = C / dAexp 
-oo « J —oc 



iXE + log (z/ + Xe p ) 



We replace the summation in the exponential by integration 



-iXE- 



gv 



W m . c . e .{v) =C j dAexp 
and expand the logarithm in the integral into the Taylor series: 

oo 

logWp [y + Xe p ) = logWp (0) + ^ 

This yields 

W m . c . e . {v)=C / dX exp -iXE + ]T - ^ 

.7 = 1 J ' n ' 



d 3 p\ogw p (v + Xe p ) 



« J ^^(^ + Aep) J 



i=i 



d^prrijip) [v + Xe p ) J 



(52) 



(53) 



(54) 



(55) 
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The first term of the sum in the exponential can be rewritten as 

. V 



(2tt) 2 



d 3 p(n p ) (y + \e p ) = i (vN + XE) , 



where N is the average total particle number (|40|l in the grand canonical system, and 

V 



E 



(2tt) 3 



dpe p (n p ) 



(56) 



(57) 



is its average total energy. The temperature of the GCE has been arbitrary so far. Now we fix it so that the average 
energy of the GCE equals the energy of the MCE: E = E. In this case the term —iXE is annihilated by iXE. Finally 
one gets 



d 3 pmj(p) (v + Xe p ) J 



(58) 



Let us find the values of the function W m . c . e .(v) and its derivatives at v — as they will be needed to calculate the 
moments of the probability distribution W m . c . e .(N). From Eq. I|58|l it follows 



/+oo 
dAexp 
-OO 

It is convenient to introduce the notation 

V 



E 



(2tt) 3 



VX 3 V 



d 3 pm 2 (p)el 



d 3 pmj{p)t 



(59) 



(60) 



After replacing the integration variable, A — x/a, the integral l|59|l takes the form 



C f +OQ ( x' 

W m .c.e, (0) = — / dxexpl — 2~)''' s]) 



j 



21 ~j 



where 



It is easy to see that 



Klj 



1 V 

a<> (2?r) 



Klj 



d 3 p mi{p) e p 



V i-i/2 



(61) 



(62) 



(63) 



i.e. the coefficients Kjj, j > 3 ,in the second exponent of (|f> 1|> become small at V — > oo. We expand the second 
exponential function in Eq. I|61|l and perform the integration: 



Wm.c.e.(0) 



5 

24 



1+ ( -K4.4-— k| i3 ) +0{V~ 2 ) 



(64) 



Similarly, 



dW m .c.e.( v ) 



dv 



i/=0 



oo •j—X 
J=2 



(7-1)1 

K 2,1 K 3,3 K 3,2 5 

2iV 



x ex P 1 — 2" ' 6XP 



v -> 2 J k: 



J=3 J 



n x 3 



2N 24 K 



3,3 



o(v- 2 ) 



(65) 
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and 



d 2 W m . c . e .(v) 



dv 2 



!/=0 



C 
a 



+ OC 



dx 



E 

i=2 



x exp ^ 


x 2 \ 
-yjexp 


E- 








K 2,l 


«3,2 , ^2,1^3^3 


" TV 2 


iV 


iV 



p 2 "Hi |iv + V 

3=2 



TV CV2¥ 



(.7-2)! 



J(J-1) j.,7-1 



(i - 1) 



^2,0 
2 



AT 



"3,3 



Now the moments of the probability distribution W m , c . e . (N) can be calculated: 



(N)r 



dW m . c . e .(v) 



and 



(N 2 ) m . c ., 



iW m .c.e (0) dv 



1 d 2 W m .c.e.H 



N 



u=0 



i + ^=(K 2 , lK3i 3-K 3:2 ) + o(y 2 ) 



Wm.c.e.(0) 



= ]v 2 



i/=0 



K 2,0 — K 2,l , «2,lK3 : 3 — K3 j2 



N 



N 



+ o(v- 2 ) 



(66) 



(67) 



(68) 



As it should be, (N) m . c , e , and (iV 2 ) m . c . e , approach, respectively, the grand canonical values N and TV in the ther- 
modynamical limit. The leading finite- volume corrections decay as 1/V. Calculation of higher-order corrections is 
straightforward. 

As is seen, the correction terms containing k 2i ik 3 .3 — K 3,2 contribute both to (N) m c e ^ and {N 2 ) m c e , but cancel 
each other in the scaling variance, 



aJ. n 



(N 2 



(N) 2 m . c . e . _ ^_4, | Q x = Jrf 3 p ^/d 3 p v 2 p e 2 p - (fd 3 p v 2 p e p ) | ^ 



(JV)m.ce. 3V ' J d^p v^J d 3 p n p 

so that the result for fluctuations is indeed the same as in the previous section (|19(l . 



IV. SUMMARY 



We have proposed a new method for a microcanonical treatment of quantum gases near a thermodynamic limit. 
The method is based on the analysis of moments of the particle number distribution in the microcanonical ensem- 
ble. For particle number fluctuations in the thermodynamic limit it leads to the same results as the microscopic 
correlator method |Tfij . However, the new method is more mathematically rigorous and consistent, and it elucidates 
some subtleties.lt gives, therefore, a justification of our previous findings [23| that the scaled variance for particle 
number fluctuations in the microcanonical ensemble is different from that in the grand canonical ensemble even in the 
thermodynamic limit. Along with fluctuations, the new method allows calculating the finite- volume corrections to the 
thermodynamic quantities in the microcanonical ensemble. This can not be done within the microscopic correlator 
calculations. Our approach can be straightforwardly extended to the system of charged particles with exact charge 
conservation laws taken into account. 
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APPENDIX A: FOURIER TRANSFORM OF A PROBABILITY DISTRIBUTION 

Let us consider a probability distribution W(x). The Fourier transform of this distribution is given by: 

W(y) = J dxe lxy W(x) , (Al) 

or, for a discrete variable x, 

W(y) = e lxy W(x) . (A2) 



It is easy to check that the derivatives of the function W(y) are related to the average values of x k : 

k 



1 fid, 

7Z W(y) 



W(0) V * d V 



The central moments of the distribution W(x) can be conveniently calculated using log W (y): 



where 



n%i = 






m 2 = 


((.*- 




m 3 = 


((.*- 




rn 4 = 


((.*- 


(*» 4 > - 


etc. 







\2\2 



(A3) 



(A4) 



(A5) 
(A6) 
(A7) 
(A8) 
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